Nonlinear atom-optical delta-kicked harmonic oscillator using a Bose-Einstein 

condensate 
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We experimentally investigate the atom-optical delta-kicked harmonic oscillator for the case of 
nonlinearity due to coUisional interactions present in a Bose-Einstein condensate. A Bose con- 
densate of rubidium atoms tightly confined in a static harmonic magnetic trap is exposed to a 
one-dimensional optical standing-wave potential that is pulsed on periodically. We focus on the 
quantum anti-resonance case for which the classical periodic behavior is simple and well under- 
stood. We show that after a small number of kicks the dynamics is dominated by dephasing of 
matter wave interference due to the finite width of the condensate's initial momentum distribution. 
In addition, we demonstrate that the nonlinear mean-field interaction in a typical harmonically 
confined Bose condensate is not sufficient to give rise to chaotic behavior. 
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The delta-kicked rotor is an extensively investigated 
system in the field of classical chaos theory. During the 
last decade great progress has been achieved in under- 
standing quantum dynamics of a classically chaotic sys- 
tem using atom-optical techniques and cold atoms. From 
an experimentalpoint of view, cold atoms in optical po- 
tentials 0, 0) 0) S H provide an ideal environment to ex- 
plore quantum dynamics. To date, all experimental work 
has focused on linear atomic systems, (see, for example, 
[^X.^.E| and references therein) where the quantum dy- 
namics is stable due to the linearity of the Schrodinger 
equation. In stark contrast to the chaotic behavior of 
classical dynamics, the linear quantum systems exhibit 
anti-resonance (periodic motion), dynamical localization 
(quasi-periodic motion) or resonant dynamics [13, [l3 • 

Recently, theoretical investigations have considered 
how the nonlinearity due to many-body (coUisional) in- 
teractions in a Bose-Einstein condensate modifies the be- 
havior of the atom-optical kicked rotor system, provid- 
ing a route to chaotic dynamics. Gardiner et al. devel- 
oped a theoretical formalism to treat the one-dimensional 
nonlinear kicked harmonic oscillator (a particular mani- 
festation of the generic delta-kicked rotor) using Gross- 
Pitaevskii and Liouville-type equations to describe the 
dynamics of a Bose-Einstein condensate, and estimated 
the grow th rate in the number of non-condensate parti- 
cles Il3| • Zhang et al. investigated the generalized quan- 
tum kicked rotor by considering a periodically kicked 
Bose condensate confined in a ring potential for the case 
of quantum anti-resonance |l3l | . As opposed to the famil- 
iar periodic behavior exhibited by a corresponding linear 
system, they predicted quasi-periodic variation in energy 
for a weak interaction strength and chaotic behavior for 
strong interactions. 

In this work we investigate the nonlinear delta-kicked 
harmonic oscillator by performing experiments on Bose- 
Einstein condensates in a harmonic potential. A Bose 
condensate of rubidium atoms tightly confined in a static 



harmonic magnetic trap is exposed to a periodically 
pulsed one-dimensional optical standing-wave potential. 
Our focus is on the particular case of quantum anti- 
resonance for which the linear behavior is simple and well 
understood 14]. The finite width of the initial conden- 
sate momentum distribution is shown to have a profound 
effect on the dynamics. After a small number of kicks the 
behavior is dominated by dephasing of matter wave in- 
terference. We present numerical solutions of the Gross- 
Pitaevskii equation which match the observed behavior 
and confirm our interpretation. 

In the atom-optical kicked harmonic oscillator, the ef- 
fective Planck's constant k can be adjusted to, in a sense, 
make the system "more" or "less" quantum mechanical. 
At specific values of ^ - in particular, where S is a rational 
multiple of 2it - quite remarkable phenomena can occur 
in the form of so-called quantum resonances and anti- 
resonances HlllIllElillillll-In this work we fo- 
cus our attention on the case of the k = 2n anti-resonance 
at which the energy of a linear system exhibits simple pe- 
riodic behavior. This anti-resonance requires a particular 
initial momentum state which evolves in such a way that 
during the period of free evolution in between kicks, the 
different components of the state vector of the system 
experience a phase shift that alternates in sign from one 
momentum component to the next, so that the system 
returns identically to its initial state after every second 
kick. The underlying physics of linear atom-optical kick- 
ing at anti-resonance has already been neatly described, 
albeit in a different context [Ijl . In the short pulse (thin 
grating) limit the first kick imprints a sinusoidal phase 
profile onto the plane matter wave thereby populating a 
number of momentum states (diffraction orders), and the 
phase evolution of the n*^ state is proportional to n'^ so 
that after free evolution (between kicks) corresponding 
to half the Talbot time {Tt = h/AEr, where the recoil 
energy Er = {hk)'^/2m, k is the wave- vector and m is 
the atomic mass) the second pulse cancels the spatial 



variation induced by the first. For multiple pulses this 
process repeats so that the initial plane wave state is re- 
constructed after every second pulse. 

Bose condensate evolution in an optical standing wave, 
or lattice, has previously been well described by the 
Gross-Pitaevskii equation (GPE) (see, for example, |21j), 
and condensate behavior in a kicked harmonic potential 
can be described in this formalism using the one dimen- 
sional GPE along the direction of the kicking beams, 
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where ip(x,t) is the condensate wave function and k ~ 
Er^TTp/h is the classical stochasticity parameter (or 
kick strength) for the effective Rabi frequency il. Here 
f{t — nT) represents a square pulse, such that f{t — 
nT) = 1 for < t — nT < Tp, where Tp is the 
pulse length. The lengt h scale is the characteristic har- 
monic oscillator length y^h/2mLUt, and the temporal scale 
is the effective trapping frequency uJt along the axis of the 
kicking beams. The nonlinear strength C = (8/i/3)^/^ is 
calculated such that the one dimensional chemical po- 
tential fi is equal to the chemical potential of the three 
dimensional condensate, in the Thomas-Fermi approxi- 
mation. Optimization techniques developed by Blakie 
and Ballagh l2l| are used to calculate the condensate 
ground state and the GPE is evolved using a Runge- 
Kutta fourth-order interaction picture algorithm J23 ■ 
The energy of the system is calculated using 
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which is evaluated after each kick to make a direct com- 
parison with experiment. 

Note that in this formalism any non-condensate par- 
ticles are not accounted for. Previous theoretical pa- 
pers [13. Il3| have investigated the proliferation of non- 
condensate particles, and for our nonlinearity, kicking 
strength, and number of kicks, this is predicted to be 
negligible. Starting with a pure Bose condensate, we do 
not observe any formation of non-condensate particles. 

Bose-Einstein condensates with up to 10^ ^^Rb atoms 
are created in the F — 2, mp — 2 hyperfine state with 
no discernible thermal component. A description of the 
BEG apparatus was given previously 23] , but there have 
been some modifications. We now use injection-seeded 
diode lasers to drive the two magneto-optical traps and 



atoms are transferred continuously between the traps 
using a focused resonant laser beam. Condensates are 
formed in the static harmonic potential of a quadrupole- 
loffe-configuration trap ^4] , characterized by radial and 
axial oscillation frequencies of ujr/2-K = 164 Hz and 
ijJz/2ti = 14 Hz, respectively. A condensate, while con- 
fined in the magnetic trap, is then exposed to a pulsed 
optical standing wave generated by two counterpropagat- 
ing laser beams with parallel linear polarizations derived 
from a single beam which is detuned 1.48 GHz from the 
5'S'i/2 7 F — 2 -^ 5P3/2, F' — i transition. Each beam has 
an intensity of 1052 W/m^ and intercepts the condensate 
at an angle of 27° to the radial direction. A double-pass 
acoustic-optic modulator is used in each beam for switch- 
ing the optical potential on and off. The pulse length is 
796 ns, which is much less than the minimum classical 
oscillation period of 130 //s, so that the kicking potential 
is well described as a thin phase grating ^&|. Following 
the kicks the momentum distribution is determined from 
a time-of-flight absorption image after a free expansion 
period of 29 ms by which time the momentum compo- 
nents have separated. The energy of the atomic sample 
is determined by calculating {J p^dp)/2m then dividing 
by the total number of atoms. The kicking period is 
33.16 fis to match the condition for the quantum anti- 
resonance at k = 8ErT/huJt = 2TT (corresponding to half 
the Talbot time), where T is the pulse period in units of 
l/u!t. The beam detuning and intensity were chosen to 
give a relatively strong kicking strength while maintain- 
ing a negligible spontaneous emission rate (< 34 s^^). 
Up to twenty five kicks were delivered to the condensate 
for each experimental run. For each number of kicks, 
the energy measurement was repeated six times, and the 
mean value is plotted in Fig. ^ 

Figure^shows an experimental and theoretical plot of 
energy versus kick number. The theoretical calculation 
uses Tp = 6.9686 x lO"-*, T = 0.029 and C = 50 in corre- 
spondence with the experimental conditions. The height 
of the optical potential has been adjusted {k = 8.25) 
so that the energy after the first kick is consistent with 
the experimental value. Initially, periodic behavior is ob- 
served, but after several kicks the oscillation in the energy 
of the system damps away to an average value (to within 
our experimental uncertainty). The theoretical points in- 
dicate that this average value gradually increases, but no 
further significant oscillation is expected. This steady in- 
crease occurs because in the time between kicks, atoms 
moving in the harmonic potential gain a small amount of 
potential energy. 

The damping in the oscillation of the energy is due 
to dephasing associated with the finite width of the con- 
densate's initial momentum distribution. The initial mo- 
mentum state is not perfectly reconstructed after each 
free evolution period because different momentum com- 
ponents of the initial distribution have a slightly different 
Talbot time (associated with their slightly different phase 
evolution). This is illustrated in Fig.|21 The rate of cou- 
pling between momentum states is not uniform across 
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FIG. 1: Energy plotted versus kick number for the quantum 
anti-resonance condition, k = 2ti with mean-field interactions. 
The solid circles are the measured mean energies and the error 
bars include shot-to-shot variation and systematic uncertainty 
in the calculation of the energy. The open circles are the 
corresponding theoretical values computed using Eq. |5| and 
the solid line is to guide the eye. 
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FIG. 2: Momentum distribution versus kick number for the 
numerical simulation in Fig.Q 



the momentum distribution of the condensate. The cen- 
tral (zero momentum) region of the initial condensate 
momentum distribution couples to the higher-order mo- 
mentum states at a slower rate than the non-zero wings 
of the condensate wave function. This causes, for exam- 
ple, the development of the double-peaked structure in 
the first order diffraction components. As time evolves 
the cycling between momentum states for different com- 
ponents of the initial distribution become progressively 



out of phase. 

This process of dephasing occurs even in the absence 
of collisions. Figure |21 illustrates the results of theoreti- 
cal simulations, showing the energy dynamics with and 
without interactions. For the experimental condition of 




10 15 

Kick number 



FIG. 3: Numerical simulations of energy versus kick num- 
ber for the quantum anti-resonance condition, k = 2ii. (a) 
No coUisional interaction (C = 0), (b) coUisional interaction 
corresponding to the experiment (C = 50), and (c) strong 
coUisional interaction (C = 1000). 

C = 50 (Fig. |2Ib)), the behavior is essentially the same 
as the collisionless case (Fig. Ofa)) although dephasing 
occurs on a longer timescale. This is due to the fact that 
the initial condensate momentum width decreases with 
increasing coUisional interactions |26|. 

After a significant time one might expect a rephasing 
of the condensate wavefunction, leading to quasi-periodic 
dynamics. While we have not observed this, our numer- 
ical simulations indicate that some rephasing is possible 
for the collisionless case, but this is highly sensitive to 
noise, and we predict that rephasing will not occur in 
the nonlinear regime. 

In Fig. OIc) the nonlinear term C in our calculation is 
a factor of 20 times larger than that corresponding to our 
experiment, and the behavior is no longer dominated by 
the damped oscillations. We begin to predict what ap- 
pears to be unstable behavior similar to that predicted by 
Zhang et al. [13, and this is consistent with our values 
of kick strength and nonlinearity. Experimentally this 
is attainable with condensates of > 10^ atoms (which 
is beyond our reach), or by increasing the s-wave scat- 
tering length via a Feshbach resonance [23 ■ Although 
Feshbach resonances have not been observed in the mag- 
netically trapped states of ^^Rb, and one would there- 
fore need to use other spin states confined in an optical 
dipole trap, this is a particularly appealing method for 
controlling the strength of the nonlinearity. Another pos- 
sibility for reaching the chaotic regime is by using a much 



lower kick strength. We repeated our measurements for a 
kick strength k — 4.125 and observed similar features to 
those presented in Fig. Q with the main difference being 
a smaller amplitude of the energy oscillations. We esti- 
mate that we would have to reduce our kick strength by 
a factor of 100 to enter the chaotic regime predicted by 
Zhang et al. [13] ■ While it may seem straightforward to 
simply further reduce the intensity of the kicking beams, 
this reduces the energy of the system to the point where 
shot-to-shot variations exceed the predicted signal. For 
a kick strength lower than k « 4 the signal to noise is 
compromised and the energy of our system becomes im- 
measurable. 

In summary, we experimentally investigated the possi- 
bility of using nonlinear coUisional interactions in a typi- 
cal Bose-Einstein condensate to observe chaotic dynamics 
in the quantum-kicked harmonic oscillator system. We 
applied a pulsed, far-detuned, optical standing wave to 
a rubidium Bose condensate, and measured the system 
energy as a function of kick number for the case of the 



quantum anti-resonance condition a,t k = 2tt. We found 
that, even in the presence of nonlinear interactions, our 
system exhibits the well-known periodic behavior associ- 
ated with the linear system. Using numerical solutions to 
the Gross-Pitaevskii equation, we showed that observed 
dephasing of the oscillations is due to the finite width 
of the condensate's initial momentum distribution. This 
severely limits the possibility of observing an extended 
period of chaotic behavior in the energy of the system. 
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